Ferromagnetism and orbital order in the two-orbital 
Hubbard model 



Katsunori Kubo 

Mcix Planck Institute for Chemical Physics of Solids, 01187 Dresden, Germany 

Advanced Science Research Center, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195, 

Japan 



Abstract. We investigate spin and orbital states of the two-orbital Hubbard model on a square 
lattice by using a variational Monte Carlo method at quarter-filling, i.e., the electron number 
per site is one. As a variational wave function, we consider a Gutzwiller projected wave function 
of a mean-field type wave function for a staggered spin and/or orbital ordered state. Then, we 
evaluate expectation value of energy for the variational wave functions by using the Monte Carlo 
method and determine the ground state. In the strong Coulomb interaction region, the ground 
state is the perfect ferromagnetic state with antiferro-orbital (AF-orbital) order. By decreasing 
the interaction, we find that the disordered state becomes the ground state. Although we have 
also considered the paramagnetic state with AF-orbital order, i.e., purely orbital ordered state, 
and partial ferromagnetic states with and without AF-orbital order, they do not become the 
ground state. 



1. Introduction 

It has been recognized that orbital degree of freedom is important for magnetism in some 
materials, such as colossal magneto-resistance and complex ordered phases of manganites [HIS]) 
and exotic magnetism in /-electron systems [21 E]- The orbital degree of freedom plays an 
important role in particular for realization of ferromagnetism. In two orbital models at quarter- 
filling, i.e., the electron number per site is one, a perfect ferromagnetic state realizes with 
antiferro-orbital (AF-orbital) order in the strong coupling limit [3]. On the other hand, in the 
intermediate coupling regime, it is not obvious whether the inclusion of the orbital degree of 
freedom is enough for the stabilization of ferromagnetism. In addition, it is interesting whether 
the perfect ferromagnetic state is continuously connected to the paramagentic state through 
partial ferromagnetic states or the ground state changes abruptly from the perfect ferromagnetic 
state to the paramagnetic state at a certain interaction strength. 

In this study, in order to unveil the ground state property in a system with the orbital degree 
of freedom in weak coupling to strong coupling cases, we study spin and orbital states of the two- 
orbital Hubbard model by using a variational Monte Carlo method at quarter-filling. A similar 
model has been studied by the variational Monte Carlo method for states without orbital order, 
and it has been found that staggered orbital correlation enhances in a ferromagnetic state in a 
strong coupling region [5]. In the strong Coulomb interaction limit, the ground state should be 
the ferromagnetic state with AF-orbital order, and in the weak coupling case, the ground state 



is the disordered state. In this study, in addition to these two states, we consider a purely AF- 
orbital ordered state without magnetic order, a purely magnetic ordered state without orbital 
order, and partial ferromagnetic state with and without AF-orbital order. 



2. Model and Method 

The two-orbital Hubbard model is given by the following Hamiltonian: 
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where Cira is the annihilation operator of the electron at site i with orbital r (= 1 or 2) and spin 
a (=1 or [), CkT-cr is the Fourier transform of it, riira = cl^^Cira, and riir = J2a ^ircr- We consider 
the nearest neighbor hopping and the dispersion is given by efc = —2t{coskx + cos ky). We set 
the lattice constant unity. The coupling constants U, U' , J, and J' denote the intra-orbital 
Coulomb, inter-orbital Coulomb, exchange, and pair-hopping interactions, respectively. 

We evaluate expectation value of energy for variational wave functions by using the Monte 
Carlo method. We consider a Gutzwiller-projected wave function of a mean-field type wave 
function as a variational wave function O EJ Ej . In this study, we consider staggered order of 
spin and/or orbital states for the mean- field type wave function. The variational wave function 
is given by 

I^) = ^gI^), (2) 

where Pq is the Gutzwiller projection operator defined in Ref [7]. The mean-field type wave 
function |$) is given by 
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where |0) is the vacuum. The quasiparticles occupy Nf, states for each spin a from the lowest 
quasiparticle energy state, where Nfj is the number of electrons with spin a. The number 
of parameters in the Gutzwiller projection can be reduced from sixteen to ten by considering 
symmetry [7], and we can further reduce the number of the parameters to seven when we fix 
Nf,. The energy of the quasiparticle in the ordered state is given by 
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The creation operators of quasiparticles are given by 
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where Q = (vr, tt) is the ordering vector considered in this study and 
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The quasiparticle gap in the ordered state is given by 

Ara = Ac + As{6^^ - 6^l) + Ao{6rl - 6^2) + ^so{Sa^ " S^l){6rl " Sr2) , (9) 

where Ac, Ag, Ao, and Ago denote the gap for charge, spin, orbital, and spin-orbital ordered 
states, respectively. We evaluate the expectation value of energy for the variational wave 
function by using the Monte Carlo method, and optimize these gap parameters and the 
Gutzwiller parameters to minimize energy. We can also evaluate energy by fixing some 
parameters, for example, we set all the gap parameters zero for the disordered phase. In 
addition to states without magnetization, we also consider states with finite magnetization 
m = {N^-N^)/iN^+N^). 

3. Result 

Figure [T] shows energy E per site as functions of U of several states for J = 0, J = O.IU, 
J = 0.2U, and J = 0.3C/ on a 12 x 12 lattice. We have used the relations U = U' + J + J' 
and J = J' . Note that at J = 0, the orbital space is equivalent to the spin space and we can 
exchange them, for example, the ferromagnetic AF-orbital state is equivalent to antiferro-spin 
ferro-orbital state. At J 7^ 0, the orbital space is not equivalent to the spin space, and the 
ground state is uniquely determined except for trivial degeneracy, e.g., the degeneracy due to 




Figure 1. Energy as functions of the Coulomb interaction obtained with several variational 
functions (a) for J = 0, (b) for J = O.IU, (c) for J = 0.2C/, and (d) for J = 0.3U. Dash-dotted 
line, dotted line, dashed line, and solid line represent energy of Hartree-Fock wave functions 
for m = 1 para-orbital, m = 1 AF-orbital, m = 0.5 para-orbital, and m = para-orbital, 
respectively. 



the rotational symmetry in the spin space. The hnes are the Hartree-Fock energy, i.e., without 
the Gutzwiher projection. 

We find that the paramagnetic (m = 0) orbital-disordered (para-orbital) state has the lowest 
energy for U < Uc, where Uc — Wt-15t. For U > Uc, the perfect ferromagnetic (m = 1) 
AF-orbital state becomes the ground state. When the Coulomb interaction becomes large, 
electrons tend to occupy the same spin state to avoid the large energy loss by the Coulomb 
interaction even though the kinetic energy becomes large. In the perfect ferromagnetic state, 
the model reduces to the single-orbital Hubbard model at half-filling with Coulomb interaction 
UeS = U' — J if we regard the spin of the single-orbital Hubbard model as orbital. Thus, in 
the perfect ferromagnetic state, orbital- AF order occurs to reduce energy. The obtained critical 
value Uc from the disordered state to the ferromagnetic state is not so large. If we do not 
consider the orbital order, Uc becomes larger. The critical value Uc increases as J is increased. 
By increasing J, the effective interaction C/eff between different orbitals becomes weak, and a 
larger value of the Coulomb interaction is necessary for the appearance of the ferromagnetic 
state. 

Although we have also calculated energy for the purely AF-orbital state with m = and 
for partial ferromagnetic states, e.g., m = 0.5 shown in figure [H with and without AF-orbital 
order, energy for them does not become the lowest. This fact indicates that the ferromagnetism 
and the AF-orbital order stabilize each other and occur simultaneously in the ground state at 
quarter-filling. Note that the absence of a partial ferromagnetic state at quarter-filling has also 
been found in one dimension [8] and in infinite dimensions [9j. 

4. Summary 

We have studied ground state of the two-orbital Hubbard model on a square lattice at quarter- 
filling within the variational wave functions. For the variational wave functions, we have 
considered Gutzwiller projected functions of spin and/or orbital ordered states. We have found 
that in the strong coupling region, the ground state is the perfect ferromagnetic state with 
AF-orbital order. In the weak coupling region, the ground state is the disordered state. We 
have also found that, at quarter-filling, other states do not become the ground state even in the 
intermediate coupling region. 
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